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ABSTRACT 

Helical magnetic background fields with adjustable pitch angle are imposed on a conducting 
fluid in a differentially rotating cylindrical container. The small-scale kinetic and current he- 
licities are calculated for various field geometries, and shown to have the opposite sign as the 
helicity of the large-scale field. These helicities and also the corresponding a-effect scale with 
the current helicity of the background field. The a-tensor is highly anisotropic as the compo- 
nents and azz have opposite signs. The amplitudes of the azimuthal a-effect computed 
with the cylindrical 3D MHD code are so small that the operation of an afi dynamo on the 
basis of the current-driven, kink-type instabilities of toroidal fields is highly questionable. In 
any case the low value of the a-effect would lead to very long growth times of a dynamo in 
the radiation zone of the Sun and early-type stars of the order of mega-years. 
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1 INTRODUCTION 



No hydromagnetic dynamo can exist driven only by differential ro- 
tation ^Elsasser|1946[ l, but it is known that such dynamos can exist 
if the turbulence is helical in the sense that its kinetic helicity 



"Hkin = {u ■ curl u) 



(1) 



(b • curl b) 



(2) 



and/or its current helicity 

'i-curr 

Mop 

do not vanish. Here u and b are the fluctuating parts of the flow U 
and magnetic field B. This condition of non- vanishing helicity is 
clearly fulfilled if the turbulence is rotating and stratified. In such 
turbulence a pseudo-scalar exists which allows the pseudo-scalars 
l[TJ and l|2j to take finite values. The same is true for linear shear 
flows where the stratified turbulence in the presence of the shear 
also can form a kinetic helicity (see Rudiger & Kitchatinov 2006). 
The simplest pseudo-scalar is the scalar product g ■ fJ with g as the 
gradient vector of the turbulence (or the fluid density) and fl the 
rotation vector. In spheres the gradient vector g is mainly radial so 
the pseudo-scalar g ■ fl has opposite signs in the two hemispheres, 
and vanishes at the equator. Because of the close relationship of 
the helicity to the a-effect in the mean-field electrodynamics of 
turbulent media. 



(the dots represent higher derivatives of Bo) the above-mentioned 
sign rules are also the sign rules of the a-effect i.e. a oc g ■ $1. 

This sort of a-effect only exists for inhomogeneous turbu- 
lence. In planetary cores, however, and also in laboratory experi- 
ments the only inhomogeneities result from boundary conditions, 
as the density gradients are negligible. One can show that under 
the presence of magnetic background fields other inhomogeneities 
also form pseudo-scalars and, as a consequence, lead to new mech- 
anisms for an a-effect (e.g. Gellert, Rudiger & Elstner 2008). In 
the present study we demonstrate that instabilities due to inhomo- 
geneous magnetic fields also lead to finite values of the helicities 
l[TJ and (|2}, and in accord with |3]( also to finite values of a. Indeed, 
in the presence of electric currents the simplest existing pseudo- 
scalar is Bo ■ curl Bo which does not vanish for helical field ge- 
ometries. We show that for such background fields the small-scale 
helicities obtain final values with the opposite sign as the helicity 
of the background field. 

According to the Rayleigh criterion, in the absence of MHD 
effects an ideal flow is stable against axisymmetric perturbations 
whenever the specific angular momentum increases outward 



^(i?W>0, 



(4) 



{u X b) = aBo + 



(3) 



where Q is the angular velocity, and (R, (f>, z) are cylindrical coor- 
dinates in a right-handed system. In the presence of an azimuthal 
magnetic field B^, this criterion is modified to 
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Bo, 4 
~R 



> 0, 



(5) 
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where po is the permeabiUty and p the density ( |Michael||1954[ l. 
Note also that this criterion is both necessary and sufficient for (ax- 
isymmetric) stabihty. In particular, all ideal flows can thus be desta- 
bilized, by azimuthal magnetic fields with the right profiles (steeply 
increasing outwards) and amplitudes. 

On the other hand, for nonaxisymmetric modes one has 



dR 



(RBl^) < 



(6) 



as the necessary and sufficient condition for stability of an ideal 
fluid at rest jTayler|1973[ l. Again, outwardly increasing fields are 
thus unstable. If Q is violated, the most unstable mode has az- 
imuthal wave numbers of m = ±1. 



2 EQUATIONS 

We are interested in the stability of the background field Bo ~ 
{0, B^{R), Bo), with Bo = const, and the flow no — 
{0, Rfl{R),0). The governing equations for the flow U and the 
field B are 



dU 

dt 



dB 

dt 



+ {U ■ V)C/ = --Vp + vMI + 

H — ^curl B X B, 
flop 



= curl {U X B)+ r]AB, 



(V) 



(8) 



and div U = div B = 0, where p is the pressure, u the kinematic 
viscosity and ri the magnetic diffusivity. Their ratio is the magnetic 
Prandtl number 



Pm 



(9) 



From now on we drop the subscripts from the large-scale values 
so that the total flow is f7 + u and the total field is B + b. The 
stationary background solution is 



il^an + —, B4, = aBH+—, 
/i^ R 

where an,bn, as and 6s are constants defined by 



(10) 



an 



1-7)2 ' 



as = 



Bin VifJ-B — V) 

R^. l-v' ' 



hi — flinRt 



bs — BinRin 
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1 — 

1 — flBV 



with 



Rii 



-ii-out ^^m 



PB = 



Bo 



Bin 



(11) 



(12) 



7?in and 7?out are the radii of the inner and outer cylinders, S7in and 
Jlout are their rotation rates, and Bin and Bout are the azimuthal 
magnetic fields at the inner and outer cylinders. A field of the form 
bn/R is generated by running an axial current only through the 
inner region R < Rin, whereas a field of the form asR is gener- 
ated by running a uniform axial current through the entire region 
7? < i?out, including the fluid. 

Given the z-component of the electric current, curU B — 
2a B, one finds for the current helicity of the background field 
B ■ curlB = 2asBo, which may be positive, negative, or zero. 

The inner value Bin is normalized by the uniform vertical 
field, i.e. 



/3 = 



Bin 

'b^ 



(13) 



For our standard profile ^s 



1 we have Bo ■ curl Bo 



2/3 Bg/pRin for the helicity of the background field. For fixed 
toroidal field amplitude this quantity scales as /3~^: 

Bo-curlBo = |. 1^. (14) 

op itin 

The sign of f} determines the sign of the current helicity. If the 
toroidal field is due to the interaction of a poloidal field with a dif- 
ferential rotation with negative shear then /3 is negative and vice 
versa. Interchanging ±/3 simply interchanges left and right spirals, 

m — >■ — m. 

As usual, the toroidal field amplitude is measured by the Hart- 
mann number and the global rotation by the Reynolds number, i.e. 

BinD „ riinB^ 



Re = 



(15) 



y^Popvrj 

D — Rout — Rin is used as the unit of length, iz/D as the unit 
of velocity and Bin as the unit of the azimuthal fields. Frequencies, 
including the rotation fl, are normalized with the inner rotation rate 
Slin. The Lundquist number S is defined by S = VPm Ha. The 
magnetic-diffusion frequency is = ri/D^ and then the Alfven 



frequency Qa 



S LJri is 



(16) 



popD 

Throughout the whole paper numerical values of helicities are 
given in units of Q\D. In this notation the helicity of the back- 
ground field can be written as 



—Bo ■ cm-lBo = ^nlD^ 

fJ-OP ip Rin 



(17) 



The boundary conditions associated with the perturbation 
equations are no-slip for u and perfectly conducting for b, at both 
R — Rin and R = iiout, where we fix 7?out = 27?in, i.e. f) = 0.5. 
The computational domain is periodic in z. The nonlinear MHD 
code used for the solution of Eqs. |7]l and l[8} has been described in 
detail by [Gellert, Riidiger & Fournier| ( |2007) (see also Foumier et 
al. 2004, 2005). 



3 RESULTS 

Fig. [T] shows the growth rates for a purely toroidal field (/? — oo), 
and no rotation. We see that beyond the critical Lundquist number, 
the growth rate is essentially linear, i.e. 



oc Q.P 



(18) 



where Pm = 1 has the steepest slope, and is thus more unstable 
than both Pm < 1 and Pm > 1. 

The azimuthal wavenumber of the modes shown in Fig.[T|is 
m = ±1. For a purely toroidal field, ±m, corresponding to left- 
and right-handed spirals, are degenerate, and necessarily have ex- 
actly the same growth rate curves. See also [Hollerbach, Teeluck &] 
|RudigerH2009[ >, who obtained the same effect in magnetorotational 
instabilities, and |Rudiger, Kitchatinov & Elstner|j201 la| l, who con- 
sider instabilities of toroidal fields in spheres. 

We next consider the nonlinear equilibration of these modes. 
As Fig. [2] shows, even though m — ±1 are degenerate, the equi- 
librated solutions do not consist of equal mixtures of both modes. 
Instead, either the left or the right mode wins out, and completely 
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Figure 1. Growth rate curves for various magnetic Prandtl numbers, for 
stationary cylinders and a purely toroidal field, with fig = 1. The growth 
rate and the magnetic field are normalized with the magnetic-dissipation 
frequency oj^ . 
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Figure 2. The left and right spirals that can be excited by a purely toroidal 
basic state. Apart from having the opposite handedness, the two modes are 
exactly equivalent; their kinetic helicity is ±6.0 ■ 10~^ and their current 
helicity is ±3.5 ■ 10^"^ (both in units of fl^D), where the left spiral is pos- 
itive and the right spiral negative. These modes do not drift in the azimuthal 
direction. The parameter values are /3 = oo, fig = !■ Re = 0, Ha = 200, 
and Pm = 1. 



suppresses the other. Which mode one obtains depends on the pre- 
cise initial conditions. If these already favor one mode, then (not 
surprisingly) that one wins, but if the initial condition is evenly bal- 
anced between the two modes, it is ultimately just numerical noise 
that determines which mode wins. Eventually though one mode al- 
ways wins; the solution consisting of an equal mixture of both is 
unstable. 

Spontaneous parity-breaking bifurcations of this type are well 
known in classical, non-magnetic Taylor-Couette flow (e.g. Hoff- 
mann et al. 2009, Altmeier et al. 2010 and reference therein), but 
are almost unknown in magnetohydrodynamic problems. To the 
best of our knowledge, the only other example is in the very recent 
work by |Chatterjee et aLlj2010| l. Given the importance of helicity 
in mean-field dynamics, any effect that generates helicity from an 
underlying basic state without helicity could be significant. 

We next present two series of solutions where Bo includes an 
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Figure 3. The turbulence intensities for flow (top) and field (bottom) fluc- 
tuations in units of D^Q^. Note that once /3 exceeds ~ 100 it has virtually 
no further influence. fj,B = 1, = 0.5, Pm = 1. 



axial component (/3 < oo). In contrast to Figs. [T| and [2] we also 
include a differential rotation here. The profiles of the basic state 
field and flow are fixed at fiB = 1 and /in = 0.5. Their amplitudes 
are Ha = 100 and Re = 200 for the first series, and Ha = 200, 
Re — 20 for the second. The first series is thus rotationally dom- 
inated (fl > Qa), whereas the second is magnetically dominated 
(SIa > f2)- The astrophysically relevant case is rotationally domi- 
nated, which is not the classical realization of the Tayler instability. 
We have called this instability the Azimuthal Magnetorotational In- 
stability (AMRI, see Hollerbach, Teeluck & Rudiger 2009). 

For both series of runs Fig. [3] shows the kinetic and magnetic 
turbulence intensities (u^) and (b^). For sufficiently large /3 its 
influence is very small; the axial component of Bq is then so weak 
that it has no further influence. This is not true for small /3, where 
the axial field starts to dominate. For /3 < 1 the kink-instability 
is strongly stabilized jRiidiger, Schultz & Elstner| ( |20Ilb} and the 
resulting energies of the perturbations are reduced. 

Fig. [4] shows the kinetic and current ^ helicities for the 
two series of runs. For both series, both helicities have the opposite 
sign as /3 (see also Kapyla & Brandenburg 2009) for comparison). 
In stating this result, it is important though to specify carefully the 
nature of the initial conditions used in each run. For /? — 0(1), 
the basic state has a sufficiently strong handedness that it forces the 
instabilities to have a particular parity as well, which as indicated 
turns out to be opposite to that of the basic state. If one then grad- 
ually increases /?, each time using the previous solution as the new 
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Figure 4. The kinetic and cuirent lielicities of the nonaxisymmetric per- 
turbations as functions of fi. Top: for f2 > SIa (Ha = 100, Re = 200). 
Bottom: for f^A > (Ha = 200, Re = 20). The dash-dotted lines indi- 
cate the limits ±6 •10"'' of the kinetic helicity of the left and right modes 
in Fig.[2] 



initial condition, this parity of the instabilities is preserved all way 
to /3 — >■ oo, where the basic state no longer has a handedness, and 
both left and right instabilities could exist equally well, as in Fig. 

El 

That is, by the time one reaches /3 = 500, say, the basic state 
makes sufficiently little distinction between left and right modes 
that both could exist, but because of the way we have reached 
(3 = 500, we consistently obtain the right mode. However, sup- 
pose one does the following experiment now: Take the right mode 
at /3 = 500, swap its parity to be left, and use that as a new initial 
condition for a series of runs where /? is now gradually reduced. 
Eventually there comes a point where the basic state's handedness 
is sufficiently great that it no longer allows the instability to have 
the 'wrong' parity, and the solution reverts back to the right mode. 
This feature that both left and right modes are allowed for suffi- 
ciently large /3 (where the degeneracy between the two modes is 
only weakly broken) but not for smaller /3 (where the degeneracy is 
strongly broken) is in many ways analogous to an imperfect pitch- 
fork bifurcation. 




Figure 5. The a-effect in the slowly rotating case, SI < Qa. Re = 20, 
Ha = 200. The top shows a^^, the bottom Ozz. The other parameters are 
/3 = 3,Ain = 0.5, MS = 1, Pm = 1. 



4 ALPHA-EFFECT AND DYNAMO THEORY 

We have also calculated the a-effect in ([3]l with the same averaging 
procedure over the azimuth. Because of the complex structure of 
the background field it is even possible to determine parts of the 
tensorial structure of the a-tensor. In particular we are interested in 
the signs and amplitudes of the a-effect in both azimuthal and axial 
directions. According to the general rule that the azimuthal a-effect 
is anticorrelated with the (kinetic) helicity we expect the azimuthal 
a-effect to be positive for /3 > 0. The expected sign of the axial a- 
effect is not clear. There are theories and simulations leading to ol^^ 
and a.zz with opposite signs (see Rudiger & Hollerbach 2004 for 
an overview). We should not be surprised to find a similar behavior 
in the present simulations. It also means that any dynamo with very 
weak differential rotation cannot be treated with a scalar a-effect. 

Figs. |5] and [6] give the results for slow and rapid rotation. On 
the basis of Eq. ol the dimensionless a-effect in the form 



Cct — 



aD 
V 



(19) 



is plotted for the components a^,^ and a^z- In both cases, these 
two components have opposite signs, with a^^ > and a^z < 
almost everywhere in the meridional plane. This anti-correlation 
between the two components is also strongest in the center of the 
gap, and weakest near the boundaries. It is therefore not caused by 
the boundaries. 
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Figure 6. The a-effect in the rapidly rotating case, D, > H a > Re = 200, 
Ha = 100. The top shows a^^, the bottom cizz - The other parameters are 
/3 = 3, /in = 0.5, MS = 1, Pm = 1. 



That Figs. [3] and [6] show such similar results is surprising, and 
is one of the basic results of this paper. The influence of rotation 
on a is evidently rather weak. The fact that - contrary to previous 
results for rotating convection - a^^ is actually smaller for rapid 
rotation than for slow rotation illustrates just how different these 
magnetic-induced helicities are from some of the previous results. 
Finally, Fig.|7]shows how the amplitudes of a^^ vary with /3, being 
roughly inversely proportional in both cases. 

To consider some possible astrophysical implications of these 
results, imagine a disk dynamo with dominant field components 
and Br. Dynamo waves of afi-type require for self-excitation 
that the product of \\9) and 



Cn 



V 



dR 



(20) 



exceeds a critical value of order unity, i.e. CaCn ^ 1. The ratio 
of the amplitudes of the field components B^ and Br follows the 
simple rule 
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Figure 7. The dimensionless dynamo number Ca of the azimuthal a-effect 
for Q < and Q > Qj^. For higher /3 (smaller B^) the o-effect de- 
creases like C/P with C ~ 0.05. /mq = 0.5, ^J.B = 1. Pm = 1- 



We know from Fig.|7]that Ca ^ C//3 with C < 1, so that p2| 
gives, at least as an order-of-magnitude estimate, the condition 



C > 



\Br 
\B. 



(23) 



for self-excitation of a dynamo with differential rotation and 
current-driven a-effect. For disk dynamos Br dominates B^, and 
for spherical dynamos Br is comparable to B^ . In both cases the 
condition for self-excitation becomes C > 1, which cannot be ful- 
filled according to Fig. [t] which suggests instead that C £ 0.05. 
The a-effect due to the current helicity of the background field ap- 
pears as much too small to allow the operation of an ail-dynamo. 

Another argument concerns the growth rate of such a dynamo 
(if it exists at all) in relation to the very long magnetic diffusion 
times in radiative zones. Assume that for self-excitation CaCn > 
1, then the growth rate cjgr is given by 



D2 



C 



\Br\ 



(24) 



Hence, only for C > 1 the growth time of the dynamo would be 
shorter than the magnetic diffusion time /rj, which is known to 
be of order Gyr for the radiative interior of stars. 

One can also argue as follows. The relation \24\ also reads 

ojgr — V aQ' (25) 

independent of the magnetic diffusivity. On the other hand, for 
given Ca {25} states 



(26) 



which for the computed value Ca — 0.01 taken from Fig. |7] and 
ri ~ 500 cm^/s for the solar core leads to values of order 10~^^ 
s~^, i.e. to growth times of order 10 Myr As it is typical for afl- 
dynamos their growth times are only slightly shorter than the basic 
magnetic decay time. 



so that dynamo excitation requires 

\B^\ 



\Bi 



\Ca > 1. 



(21) 



(22) 



5 SUMMARY 

We have shown that the current-driven instability of helical large- 
scale fields does produce small-scale helicity (kinetic plus current 
helicity) and even a-effects, but the resulting numerical values 
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seem to be too small for the operation of large-scale dynamos in 
radiative zones of early-type stars. 
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